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Preface 


With  the  advent  of  sophisticated  ballistic  missile 
early  learning  and  tracking  systems,  the  fact  that  a missile *s 
trajectory  can  be  measured  and  predicted  allows  for  a rela- 
tively easy  high  altitude  intercept.  It  has  been  suggested 
that  a modification  from  the  original  trajectory  would 
improve  vehicle  survivability.  My  goal  was  to  investigate 
this  type  of  maneuver  and  to  limit  the  numerous  available 
modifications  by  an  application  of  optimization  for  parameters 
of  interest  in  the  ballistic  missile  problem. 

I would  like  to  thank  my  advisor.  Major  Gerald  K. 
Anderson  for  his  helpful  suggestions  and  constructive 
comments  for  this  study.  I would  also  like  to  acknowledge 
Captain  Richard  K.  Potter  for  introducing  me  to  the  numerical 
methods  of  ballistic  trajectory  calculations.  Finally, 

I would  lik:  to  recognize  my  wife,  Susan,  for  her  patience, 
understanding,  and  encouragement  during  this  study. 
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Notation 


a seiiii-EaJor  axis 

DU  earth  canonical  distance  unit  = 20925740  ft 

DU/TU  earth  canonical  -elocity  unit  = 25936,28  ft/sec 

e eccentricity 

E eccentric  anocaly 

$ specific  sechanical  energy 

? hit  eqsiation 

h specif  igular  ffiomentua 

h^Q  biimou  t altitude 

h^g  reentry  altitude 

J cost  function 

J augnenced  cost  function 

P direction  to  p*erigee  in  the  perifocal  frase 


non -4  i aen s 1 ona 1 


rexocitv  laraaeter 


perifocal  direction  p-erpendicular  to  P 


radius  rector  i^roiitude 


radius  rector  s^snitu 


;he  transfer  point 


earth  concnical  tine  unit  = 806,8136  ss 


•eaction  tiae 


I 


:iae  of  flisrht  between 


:ut  ana  anojsee 


rei oclty 


cost  function  weighting:  factor 


inertial 


Symbol 

X 


Inertial  X direction  position  component 

Y inertial  Y direction  position  component 

Z inertial  Z direction  position  component 

AV  ts^nsfer  velocity  Impulse 

A Lagrange  multiplier 

3 2 

gravitational  parameter  1 DU  /TU 

V true  anomaly 

true  anomaly  at  transfer 
flight  path  angle  at  burnout 
U/  free  flight  range  angle 


Superscripts 
- vector 

T transposed  vector 


Subscripts 

a apogee 

bo  burnout  point 

f vehicle  at  final  time 

ff  free  flight 

m modified  orbit 

mod  trajectory  modification  point 

n nominal  orbit 

p perifocal  P direction 

p*  modified  orbit  component,  written  in  nominal 

perifocal  P direction 

q perifocal  Q direction 

vill 


S 


Subscripts 


modified  orbit  component,  written  in  nominal 
perifocal  Q direction 

reentry  point 

required  value 

target  at  final  time 

inertial  X direction  component 

inertial  Y direction  component 

inertial  Z direction  component 

vector  component 

vector  component 

vector  component 
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Abstract 


^he  problem  of  finding  optimal  mid-course  modifications 
of  ballistic  missile  trajectories  is  investigated,  A single 
velocity  Impulse  is  applied  at  the  point  of  transfer  from 
the  original  orbit  to  the  modified  trajectory.  This  limits 
the  modified  trajectory  to  one  which  intersects  the  original 
trajectory.  Optimization  of  a weighted  function  of  transfer 
velocity  impulse  and  che  time  to  impact  after  the  modification 
occurs  is  accomplished.  The  function  can  be  weighted  to 
accommodate  trade-offs  between  both  components. 

The  study  is  divided  into  two  parts.  Elliptical  orbits 
for  a non-rotating,  atmosphere  free  earth  are  investigated, 
and  then  trajectories  for  an  oblate,  rotating  earth  with 
atmospheric  reentry  are  examined.  For  the  elliptical  orbit 
cases,  modified  and  nominal  orbits  are  coplanar.  Pre-apogee, 
apogee,  and  post-apogee  transfers  from  the  nominal  trajectory 
to  a modified  orbit  at  apogee  are  coi.sidered.  It  is  found 
that  a pre-apogee  transfer  from  a near  circular  nominal 
trajectory  is  advantageous  for  the  defined  problem. 

A method  of  computing  an  optimal  transfer  from  a single 
point  on  a lofted  nominal  trajectory  for  nonplanar  nominal 
and  modified  trajectories  is  then  presented.  An  algorithm 
is  derived  where  a nominal  trajectory  to  a pseudo  target  is 
calculated,  then  transfer  to  a modified  trajectory  which 
imi»cts  the  real  target  is  computed.  An  example  problem  is 
shown  and  the  region  of  a minimum  cost  transfer  is  found. 
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OPTIHilL  MID-COURSE  MODIFICATIONS 
OP  BALLISTIC  MISSILE  TRAJECTORIES 


I.  INTRODUCTIOK 


Background 

A long-range  ballistic  Filssile  trajectory  can  be 
divided  into  three  phases,  powered  flight,  free  flight, 
and  reentry.  After  the  initial  powered  flight  ceases,  the 
entire  trajectory  of  the  vehicle.  Including  impact  point  can 
be  calculated  through  use  of  a series  of  radar  measurem rats. 

Green  (Ref  4i6)  comments  on  the  anti -ballistic  missil^e 
Intercept  problem  during  the  free  flight  phase i Once  the 
launch  has  been  detected,  early  warning  and  tracking  infor- 
mation can  be  supplied  to  the  free  flight  tracking  system  of 
a defender.  This  mid-course  tracking  system  needf  to  search 
only  a small  space-time  zona  to  acquire  the  vehicle  for  trajec- 
tory prediction.  The  long  free  flight  phase  is  then  a 
disadvantage  to  the  attacking  vehicle  because  l.utercept  is 
easier. 

To  counteract  this  prospect  of  mission  failure,  it  may 
be  advantageous  to  make  the  reentry  vehicle  c’range  course 
and  target  at  some  time  during  free  flight.  However,  other 
factors,  such  as  accuracy,  additional  mass  required,  cost, 
and  reliability  must  be  considered.  If  an  acceptible  mid-course 


modification  system  is  employed,  chances  of  vehicle 
survival  are  increased  because  the  missile  Is  not  committed 
to  a single  trajectory  after  powered  flight  ends. 

Bamaby  (Ref  1j16)  suggests  that  this  mid-course 
modification  is  a viable  penetration  aid  for  use  against 
high  altitude  intercept  of  anti-ballistlc  missiles.  Addi- 
tionally, if  this  maneuver  could  be  masked,  for  example,  by 
releasing  chaff  at  the  point  of  trajectory  change,  all 
defensive  radar  contact  might  be  lost. 

In  any  case,  the  defender  has  less  time  for  trajectory 
calculation  if  the  vehicle  is  again  acquired  by  radar.  This 
reduction  in  defensive  reaction  time  means  that  a less 
accurate  measurement  of  the  trajectory  is  made  and  the  high 
altitude  ABM  threat  may  be  circumvented. 


Problem  Statement 

The  purpose  of  this  study  is  to  investigate  mid-course 
modifications  of  long-range  ballistic  missile  trajectories. 
Optimal  orbit  changes  are  found,  where  the  optimality  of 
the  new  orbit  is  defined  in  terms  of  the  total  velocity 
impulse  required  to  change  the  trajectory,  and  the  warning 
time  given  to  the  defender  after  the  ti'ajectory  is  modified. 

The  ballistic  ti^Jectory  modification,  as  a specialized 
case  of  the  general  orbital  transfer  problem,  can  be  accom- 
plished through  any  number  of  thrusts.  This  investigation 
limits  the  transfer  to  a single  impulsive  thrusting  manuever, 
which  in  turn,  requires  that  the  original  and  modified 


trajectories  Intersect. 
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Outline 


I 


The  mid-course  modification  problem  is  developed  in 
two  stages.  First,  the  two-dimensional  transfer,  where 
the  nominal  and  modified  orbits  are  ccplanar,  is  examined. 

In  Chapter  II  this  type  of  transfer  is  defined.  The 
ncninal  t.-:ajectory  parameters  are  developed  in  Chapter  III, 
Chapters  IV,  V,  and  VI  present  three  types  of  coplanar 
tra.3ectory  modifications,  along  with  the  optimum  results 
for  each. 

Chapter  VII  poses  the  general  thr-ie-dimensional  or 
nonplanar  trajectory  change  problem.  Trajectory  modifications 
and  optimization  for  this  case  are  presented  in  Chapter  VIII, 
Tnls  chapter  also  discusses  results  of  the  nonplanar  modifi- 
cation case.  Finally,  the  conclusions  and  recommendations 
resulting  from  this  study  are  presented  in  Chapter  IX, 
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II,  PBOBLEM  DEFINITION  - COPIANAP  TRAJECTOBISS 


In  this  chapter,  the  l^sic  foundation  for  the  ooplanar 
transfer  problem  is  established.  Problem  definition  Inclines 
the  assumptions  made  to  analyre  the  transfer  problem,  the 
type  of  nominal  traJectoiT  used  in  this  part  of  the  study, 
the  kinds  of  orbital  transfers  considered,  and  a discussion 
of  optimization  pertaining  to  the  trajectory  modification 
problem. 


Assumptions 

For  coplanar  nomlral  and  modified  trajectories,  a 
spherical,  non*rotating  earth  ifith  no  atmosphere  is  assumed. 
The  equations  of  a oonlc,  specifically,  the  eqtiation  of  an 
ellipse,  then  describes  the  ballistic  trajectory.  This 
investigation  did  net  Include  transfer  to  parabolic  and 
hyperbolic  trajectories  since  large  velocity  impulses  are 
ne^ed  to  establish  these  types  of  orbits, 

Bsu:^b  canonical  units  are  used  in  the  analytic 
expressions  describing  distance,  time,  and  velocity.  These 
units  are  defined  by  Bate  (Ref  2i429): 

distance  unit  1 DD  = 20925740  ft 

time  unit  1 TO  « 806,8136  sec 

speed  unit  1 KF/TO  > 25936.28  ft/seo 

gi^vltational  j^rameter  1 DoVto^  * 1,407654x10^^  ft^/sec^. 

Target  location  and  missile  altitude  at  burnout  are 
two  basic  parameters  in  the  caloulatiaas.  The  target  rsaig^ 


is  aBStan^  to  be  a typical  ICEM  range  of  6000  n,  mi,  Bnmont 
is  assumed  to  occur  av  ,05  DU  altitude. 

These  initial  assumptions  can  be  used  to  solve  the 
transfer  problem  in  closed  form,  If  the  correct  nominal 
trajectory  i3  employed, 

nominal  Trajectory  Definition 

There  are  three  different  types  of  ballistic  missile 
trajectories,  the  low,  the  high  or  lofted,  and  the  maximum 
range  trajectory.  Both  the  lofted  and  the  low  trajectory 
can  reach  a target  with  the  same  initial  conditions  of  veloc- 
ity and  altitude  at  burnout.  The  distinguishing  factor  is 
the  direction  of  the  velocity  vector  at  burnout. 

The  low  ballistic  trajectory  is  characterized  by  its 
smaller  apogee.  This  type  of  orbit  has  the  advantage  of 
r^ichlng  a target  quickly,  since  less  total  distance  along 
the  trajectory  is  covered*  however,  this  type  of  trajectory 
is  not  generally  useful  for  ballistic  missiles  since  the 
vehicle  nay  encounter  the  a^SlOsphe^e  too  quickly.  Even  though 
this  section  does  not  consider  atmospheric  reentry,  the  low 
trajectory  is  not  used  since  it  is  not  a realistic  situation 
for  lallistlo  missiles. 

The  lofted  trajectory  is  produced  by  teving  a higher 
elevation  of  the  velocity  vector  at  burnout.  This  path 
las  the  advantage  of  being  less  susceptible  to  initial  errors 
at  burnout.  Since  this  trajectory  has  a higher  apogee,  more 
distance  is  covered  than  for  a low  trajectory  to  the  same  target. 
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As  would  be  expected,  the  time  of  flight  is  longer  for  this 
case.  Because  the  final  accuracy  is  better,  this  type  of 
ti^jectory  is  routinely  employed  for  ballistic  missiles.  The 
lofted  nominal  trajectory  is  used  for  the  second  part  of  this 
study  and  is  discussed  further  in  Chapter  VII, 

The  third  type  of  trajectory,  the  one  for  maximum  range, 
is  the  lowest  energy  orbit  needed  to  reach  the  target.  For 
given  initial  burnout  conditions  there  is  usually  a high  and 
low  trajectory  to  a given  target  range,  but  there  is  one  target 
range  where  only  one  path  to  the  target  exists.  This  is  the 
maximum  range  trajectory. 

Maximum  range  trajectories  exist  only  for  target  ranges 
below  180°,  For  a range  of  180°,  the  maximum  range  trajectory 
is  a circular  orbit  if  burnout  occurs  above  the  earth's  sur- 
face, The  vehicle  will  not  Impact  on  the  earth's  sxirface  for 
this  range.  If  burnout  occurs  at  sea  level,  this  trajectory 
is  a path  along  the  surface  of  the  earth  and  cannot  be  used. 

If  the  maximum  range  trajectory  is  assumed  symmetrical, 
i,e,,  burnout  and  reentry  points  are  at  equal  altitudes,  then 
the  free  flight  portion  of  the  orbit  can  be  completely  spec- 
ified by  the  free  flight  range,  or  the  range  angle  subtend«l 
at  the  ^rth's  center,  and  the  burnout  altitude. 

Because  this  maximum  range  trajectory  gives  the  most 
range  for  an  initial  velocity,  and  since  it  can  be  analyt- 
ically expressed  by  two  parameters,  it  is  selected  for  use 
as  a nominal  trajectory  in  the  coplanar  t? insf er  problem. 
Nominal  naxlmum  range  trajectory  jai^meters  are  developed 
in  Chapter  III, 
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gholoe  of  Modified  Orbits 


In  addition  to  defining  the  nominal  traJectOT7,  the 
modified  orbit  i^ch  cemneets  a transfer  point  on  the  nominal 
to  a target  Is  required.  A general  transfer  at  any  point  in 
bho  nominal  trajectory  includes  all  types  of  modified  orbits. 
Sereral,  less  complex  analytic  elliptical  cases  are  e:»minM 
in  this  study.  Transfer  to  the  apogee  of  a modified  orbit  is 
a logical  choice  since  the  Telocity  oKignltikle  at  apogee  is  the 
lowest  for  the  modified  trajectory.  The  type  of  modified 
trajectories  studied  are  as  follows t 

1,  Transfer  from  a nominal  trajectory  at  apogee  to  a 
new  orbit  idiich  also  i»B  its  apogee  at  the  transfer  point. 

2,  Pre-apogee  transfer  from  a nominal  trajectory  to 
a now  orbit  which  l»s  its  apogee  at  the  transfer  point, 

3,  Post-apogee  transfer  from  a ntadnal  trajectory 
to  a new  orbit  which  has  its  apogee  at  the  transfer  point. 

These  coplanar  transfers  are  develops!  in  Chapter  IV, 

V,  a»l  VI  respectively.  The  equations  of  an  ellipse  are  usai 
a^  all  orbital  parameters  are  found  analytically. 

Tra  jeotory  Optimization 

A number  of  modlfi^  orbits  satisfy  the  iHPOblem  of 
impacting  upon  a target  from  some  point  in  the  nominal  tra- 
jectory. In  OTder  to  relate  the  merits  of  me  trajectory  to 
t)»8e  of  another,  a cost  function  must  be  defined. 

The  factors  considered  in  the  cost  fraictlafi  formulation 
are  the  impulse  required  to  transfer  to  the  medlflGd  orbit 
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awl  tha  raaetlcm  tiae  given  to  a defender  after  this  change 
Is  Bsde,  It  Is  assuaed  that  the  booster  Is  capable  of 
lapartlng  velocity  needed  to  establish  the  noainal  trajectory 
so  this  is  not  considered  in  the  cost. 

The  iapnlse  required  to  transfer  from  the  noainal  to 
the  aodiflM  trajectory  is  directly  related  to  the  aacunt 
of*  propellant  used  by  the  isropulsion  systea.  The  transfer 
iapulse  is  defined  at  the  point  of  transfer  as  the  asgnitude 
of  the  difference  between  the  two  different  trajectory 
velocities i 

AV=|V^-VJ  (1) 

where  is  the  velocity  on  the  nodifled  orbit  at  the 
transfer  point  and  Vj^  is  the  velocity  on  the  noainal  trajec- 
tory at  the  transfer  point. 

Reaction  tiae  is  included  in  the  cost  foraulation  since 
vehicle  survivability  is  related  to  the  aaount  of  warning 
given  to  the  defender.  But  since  both  tiae  smd  A7  are 
Included  In  the  cost  foraulation,  soae  type  of  relation 
between  the  two  components  aust  be  consider^. 

To  produce  this  relationship,  a weighting  factor  on 
reaction  time  is  Introduced,  For  high  values  of  this  factor, 
the  ruction  time  contributes  aost  to  the  cost.  Low  values 
of  the  factor  allow  transfer  AV  to  dominate.  The  cost  fimctlon 
is  then  defined  as 
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/•ill, 


J = |av2*^t2 


(2) 


where  AV  is  the  transfer  lapulse 
V is  the  weightins  factor 
Tj.  is  the  reeotit^  tise. 

The  feature  of  incltiding  ¥ in  the  cost  function  allows 
a trajectory  aesigner  to  apply  the  optiaization  to  different 
■id-course  propulsion  systems.  The  effect  of  rarying  W ms 
ezaained  at^  the  walue  of  ,15  was  selected  since  this  gare 
optimum  transfers  in  a realistic  range  of  transfer  AY.  The 
trajectories  i^ch  produce  the  lowest  values  of  this  cost 
functicm  are  the  desired  soluixons  to  the  defined  irohlea. 


III. 


MAXIMUM  RANGE  NOMINAL  TRAJECTORY  CALCULATION 


Noalml  Orbit  ^raaetcrs 

If  a naxliHuia  range  trajectory  is  assumed  to  be 
synnaetrlcal,  i,*,»  burnout  and  reentry  occur  at  the  same 
altitude,  the  sagnltude  of  the  radius  vector  at  burnout, 
r^Q,  and  the  free  flight  range  angle,  describe  all  elements 
of  the  two-dimensional  trajectory.  This  nominal  trajectory 
is  shown  lii  fig,  1.  It  should  be  noted  that  is  the  angle 
* a tween  the  burnout  and  reentry  radii. 

The  flight  path  angle  which  relates  local  horizontal 
to  the  local  velocity  vector  is  defined  for  the  maximum  range 
trajectory  (Ref  2*292)  as 


<P 


be  4 


(3) 


where  ^bo  is  the  flight  path  angle 

is  tne  nominal  free  flight  range  angle. 

The  nominal  orbit  eccentricity  can  be  derived  from 


2€h 


2] 


1/2 


where  e^j  Is  the  nominal  orbit  eccentricity 

S is  the  specific  mechanical  energy  of  the  orbit 


h is  the  specific  angular  moaentun  of  the  orbit 
p is  the  earth  gravitational  parameter. 


Specific  mechanical  energy,  constant  for  the  orbit, 
Is  found  from  the  relation 


e- 


Vbo  if: 

2 


(5) 


r 

'bo 


where  Vbo  is  the  burnout  velocity. 

The  specific  angular  momentum,  also  constant  for  the 
nominal  tra  jector.; , can  be  found  at  burnout » 


f'  = ''boVboCos<fcbo 


(6) 


Equations  (5)  and  (6)  require  the  value  of 
which  can  be  found  from  the  nondlmenslonal  orbital 
parameter,  defined  by  Bate  (Ref  2t280)j 

-v/2 


^bo"  '^bo  '"bo 


(7) 


This  parameter  is  the  squared  ratio  of  the  velocity  of 
the  vehicle  to  circular  orbital  velocity  at  the  burnout 
point.  For  a oailmum  range  trajectory  (Ref  2i293) 


Q 


bo 


, 2 S\n{\/2) 

1*Sin(Y  /2) 


(8) 
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By  substituting  eqs  (5)»  (6),  and  (7)  into  eq  (4), 
the  following  expression  for  nominal  orbit  eccentricity 
results* 


n"  J +(Qbo“^^ ^bo- 


Thls  can  be  simplified  to  a function  of  and  by 
substituting  eq  (8)  into  eq  (9)» 


- j 1 - ‘^Sin(^p,/2)Cos^4>bo 

■ 1 ~[l*Sin(V2)F 


1/2 


Nominal  Time  of  Free  Flight  Calculation 


Time  of  flight  on  the  elliptical  orbit  Is  found  through 
use  of  the  Kepler  time  of  flight  caloulaticns  (Ref  2il85).  The 
true  anomaly  at  burnout  for  the  maximum  range  trajectory  is 


(11) 


where  is  the  true  anomaly  at  burnout. 

Once  knowing  t'to*  eccentric  anomaly,  S^q»  ®ay  be  found* 


P^--Cos-'l^n*CosUi,Q  \ 


(12) 


The  nominal  orbit  semi -major  axis,  a^^.  Is  found  through 
use  of  eqs  (5)f  (7)t  atid  (8)  substituted  Into 


C-- 


H- 


23^ 


(13) 


This  yields  the  desired  result  for  ajji 


n 


2 


(14) 


Sinoe  Sjj  and  are  knovm,  the  nominal  trajectory  free 
flight  time  can  be  found.  Because  the  orbit  Is  assumed  to  be 
symmetrical,  the  total  time  of  flight,  is  twice  the  time 

from  burnout  to  apogee,  and  by  Inspection  of  fig,  1,  the 
eccentric  anomaly  at  apogee  Is  tt.  Then 


3/2 


('TT-Ej^-^epSinEK^) 


bo 


(15) 


It  has  been  shown  that  all  necessary  pai^meters  for 
the  nominal  maximum  range  trajectory  can  be  found  as 
functions  of  the  bumout  i^dlus,  rt>o,  and  the  free  flight 
range  angle,  tj/^. 
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IV.  APOGEE  TO  APOGEE  TRAUSFEH3 


The  first  case  considered  is  the  transfer  from  a nominal 
orbit  at  apogee  to  a new  orbit  also  at  apoge.e.  This  type  of 
transfer  is  depicted  in  fig.  2,  Since  the  e.pcgee  radius  of 
the  nominal  trajectory  defines  the  transfer  point  and  the 
equations  of  an  ellipse  are  employed,  the  t3*ansfer  is  only  a 
function  of  the  desired,  final  range  angle  and  the  eccentricity 
of  the  modified  orbit. 

The  required  modified  Tange  angle  emplcyed  In  all  analytic 
examples  is  100°,  equivalent  to  6000  nautical  miles  on  a 
spherical  earth  surface.  The  nominal  orbit  used  is  defined 
in  Chapter  III, 

Modified  Orbit  Parameters 

The  semi -major  axis  of  the  modified  orbit  is  defined  by 

~ (16/ 

idiere  is  the  modified  orbit  semi-sajor  axis 
is  the  nominal  orbit  apogee  radius 
e^  is  the  modified  orbit  eccentricity. 

Through  the  use  of  the  equation  of  a conic 

p,  a(l-e^) 

1 + eCos  u 


(17) 
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any  true  anomaly,  p , can  be  found  in  terms  of  r,  a,  and  s. 

It  is  necessary  to  find  the  true  anomaly  at  the  point  where  ■ 
the  modified  orbit  impacts  the  earth  in  order  to  obtain 
the  final  range  angle  resulting  from  the  orbit  modification. 

Using  eq  (17)  and  solving  for  the  true  anomaly  at  the 
target, 


U.  = 


Cos  ^ { 1 


m 


3rv-|(l  - ) 


- 1 


(18) 


where  r^  Is  the  radius  at  the  target  location.  It  must  be 
noted  that  lies  between  wand  2w, 

For  impact  at  the  target  on  the  spherical  earth,  the 
value  of  r^.  is  1 DU,  Then 


a (1-e^)- 
m m 


(19) 


The  final  range  angle  is  calculated  once  the  values  of 
true  anomaly  at  the  burnout  and  Impact  points  are  known.  The 
range  angle  traversed  on  the  nominal  trajectory  is  from  burn- 
out to  apogee  or  w-  The  range  angle  traversed  on  the 

modified  orbit  is  the  angular  difference  between  apogee  and 
Impact,  Combining  these  two  ranges,  the  resulting  range 

angle,  is 
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(20) 


Cost  Function  Barametgrs 

The  tifo  parameters  required  to  define  a cost  for  each 
transfer  are  fotind  analytically.  Time  of  flight  calculations 
are  done  using  the  Kepler  time  of  flight  eqiiatlons  (Ref  2il86), 
Time  of  flight  on  the  nominal  trajectory  is  one-half  the 
nominal  tine  of  flight  derived  in  Chapter  IIIi 


'f12‘®n  '‘-bo*®n^'b  E[jq) 


(21) 


vhere  T^j^2  flight  from  burnout  to  nominal 

orbit  apogee  (TU), 

To  find  the  time  of  flight  from  apogee  to  Impact,  or 
the  reaction  time,  it  is  first  necessary  to  find  the  eccentric 
anomaly  at  Impact,  Ef 


Cos'’ 


/em-Cosv^  \ 


(22) 


also  lies  between  v and  Btt, 
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Now,  the  reaction  time  can  be  foimdi 


Tp=aj^^^(Et-enr|Sin  E^’-rr) 


(23) 


where  Tj.  Is  the  reaction  time. 


To  find  the  velocity  impulse  at  transfer,  the  apogee 


velocities  of  both  orbits  are  calculated  and  the  difference 


computed,  A perifocal  coordinate  system.  Illustrated  In 


fig,  2,  with  P In  the  direction  from  the  orbit  focus  to 


perigee,  and  Q perpendicular  to  P In  the  direction  of  travel 


on  the  orbit,  is  defined. 


Velocities  in  the  perifocal  frame  at  any  point  can  be 


foimd,  once  orbit  parameters  are  known  (Ref  2i73)i 


V=  — [-Sim/P  + (Cosi/*e)Q]  ^24) 

r(l^eCosi^) 


At  apogee,  v = ir  which  means  that  Vp  = 0,  Velocity  in  the 


Q direction  Is  then 


V_=- 


1-e 


(25) 


The  AV  required  at  apogee  can  then  be  calculated! 


AV=|V  -V  I 
• qm  qn‘ 


1-e 


n/2 


m 


1- 


\1/2 


n 


(26) 


Since  AV  and  T^.  are  known,  a cost  may  be  found  for  each 
nominal  trajectory  and  its  corresponding  solution. 

Method  of  Solution 

For  the  apogee  to  apogee  transfer  case,  it  is  desired  to 
find  a modified  orbit  eccentricity  which,  for  a given  nominal 
trajectory,  hits  a specified  range  angle,  The  method  of 

solution  involves  one  function  of  the  variable  ejjj.  The 
function  is  the  "hit"  equation,  Fi 


F=0='^ 

^ r 


(27) 


When  F=0,  the  modified  orbit  Impacts  at  the  required  target 
range,  Finding  the  solution  to  eq  (2?)  involves  an 
iteration  on  em.  A general  purpose  algorithm,  presented 
by  Powell  (Ref  7),  solves  a system  of  nonlinear  equations. 


This  method  is  used  for  the  iteration  on  Bjj,  and  other 
iterative  formulations  derived  in  this  investigation.  A 
further  dis.'ussion  of  Powell’s  algorithm  and  its  use  is 
found  in  Appendix  B. 
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For  a given  nominal  trajectory,  only  one  solution  of 


eq  (27)  exists  In  this  type  of  transfer. ♦ A series  of 
solutions  Is  geneiated  by  varying  and  a cost  function 
is  formed  from  the  set  of  solutions.  The  computer  program 
used  for  this  Is  listed  as  Program  1 in  Appendix  C, 

Since  the  cost  function  is  bounded  - the  lower  limit 
is  zero  and  the  upper  limit  is  within  the  i»rabolic  modified 
orbit  limitation  - a minimum  cost  exists.  For  a minimum  cost 
solution,  an  associated  nominal  trajectory  can  be  found.  This 
nominal  trajectory  is  then  used  as  the  initial  targeting 
information  for  the  optimal  apogee  to  apogee  transfer*  >o 
the  minimum  cost  solution  must  be  found  to  define  the  nomlr»l 
trajectory  necessary  for  an  optimal  apogee  to  apogee  transfer. 

For  the  problem  formulated,  an  interior  minimum  is 
exhibited.  The  minimum  cost  value  and  the  associated  nominal 
trajectory  can  be  found  by  considering  a constrained 
optimization  problem.  The  cost  function,  eq  (2),  is  augmented 
with  the  constraint  that  the  hit  equation  be  solved; 


2 2 ^ m r 


(28) 


where  A is  the  Lagrange  multiplier  associated  with  the  hit 
condition. 


♦ The  other  coplanar  transfers  considered  have  an  infinite 
number  of  solutions  for  each  nominal  trajectory. 
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will 


The  minimum  of  eq  (28)  is  found  by  solving  the  first 
order  necessary  conditions  of  a mlniBmm.  Three  variables, 
®m*  *^n»  ^ following  equations  to  find 
the  minimum  point: 


dl 

dV, 


= 0=Av 


n 


6AV 

6V, 


•WT, 


6T, 


X 


m 


n 


6Vn  6V. 


n 


(29) 


dj 


=0  = AV 


'm 


6e 


m 


X - rn 


(30) 


^ J “0  ~ ^ ^ 

6X'  ’ 


(31) 


To  solve  eqs  (29),  (30),  and  (31)#  the  partial  derivatives  of 
J with  respect  to  and  e«  are  required.  These  derivatives 
are  calculated  through  use  of  a first  difference  scheme 
(Ref  3*217),  It  was  found  that  regulating  the  initial 

perturtetions  so  that  first  differences  are  on  the  order  of 

-8 

10  allows  an  accurate  minimus  point  to  be  calculated. 

This  snail  firsu  difference  is  well  within  the  15  significant 
digit  capability  of  the  computer  used. 

The  first  order  necessary  conditions  of  a minimum  were 
progs^mmed  and  solved  for  this  case.  A listing  of  the 
computer  program  is  foxmd  as  Progias  2 in  Appendix  C. 
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Results  and  Analysis 

For  a target  range  of  100°,  solutions  for  transfer 
from  noalnal  range  angles  between  70°  and  180°  are  found, 
using  PrograjB  1,  An  interior  ainlHUS!  In  the  cost  function 
is  c^diibited,  and  the  alnisraa  point  is  found  by  employing 
Program  2. 

Figure  3 shows  the  resulting  cost  fianction  for  a 
weighting  factor  of  .15  . Figures  4 and  5 depict  the 
ruction  time  aiKi  velocity  which  prcxiuce  the  cost  function. 

A siniHua  cost  transfer  for  the  apogee  to  apogee 
transfer  case  occurs  for  a nominal  range  angle  of  145.85°. 
Since  the  altitude  at  burnout  is  assumed  to  be  .05  DU,  this 
nominal  range  angle  is  sufficient  to  define  al?.  initial 
burnout  parameters  of  the  nominal  trajectory.  The  optimal 
apogee  to  apogee  trajectory  sMifl cation  is  then  produced 
from  this  initial  trajectory. 

The  behavior  of  the  cost  function  for  this  case  can 
be  explained  by  analyzing  figs.  4 and  5.  As  nominal  range 
incr^ses,  a larger  mt^ifiea  orbit  eccentricity  is  required 
in  order  to  hit  the  target  range.  The  higher  eccentric 
orbits  cause  a r«iuction  in  reaction  time,  as  seen  in 
fig.  4.  This  effect  is  shown  in  eq  (23). 

Because  transfer  to  a higher  eccentric  orbit  is 
required  for  long  nominal  range  angles,  the  velocity  needed 
to  transfer  to  the  modified  orbit  incr^ses.  This  is  shown 
in  fig.  5. 
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PRE-APOGEE  0?O  APOGEE  TRANSFERS 


The  second  case  considered  is  the  pre-apogee  transfer 
from  a nominal  orbit  to  apogee  on  the  new  orbit.  A limiting 
case  of  this  particular  transfer  is  the  apogee  to  apogee 
transfer  described  in  Chapter  IV.  The  fact  that  transfer  may- 
be accomplished  prior  to  nominal  apogee  introduces  an  added 
degree  of  freedom  from  the  first  case.  Three  parameters 
are  now  required  to  define  the  modified  orbit.  They  are  the 
modified  orbit  eccentricity,  Cjjj,  the  radius  at  transfer,  r2» 
and  the  desired  final  range  angle,  The  pre-apogee  transfer 

is  shown  in  fig,  6, 

This  transfer  problem  Is  very  similar  to  the  first  case 
considered.  The  only  differences  arise  in  the  definition  of 
the  transfer  point,  calculation  of  the  modified  range  angle, 
and  finding  the  transfer  AV,  The  nominal  orbit  is  the  maximum 
range  trajectory  described  in  Chapter  III, 

Modified  Orbit  Parameters 

The  orbital  parameters  for  this  case  are  those  found 
in  Chapter  IV,  with  the  following  exceptions.  Definition 
of  the  transfer  point  on  the  nominal  trajectory  involves  the 
calculation  of  a true  anomaly  at  transfer,  ^2*  Because  the 
radius  at  transfer,  r2»  is  variable,  the  true  anomaly  at 
transfer  can  be  found  from  eq  ( i? ) s 


2? 


•Apogftft  to  Apogee  Transfer 


(32) 


i^2  = Cos“^ 


) 


Since  this  case  considers  pre-apogee  transfers, 
between  0 and  tt.  The  variable  r^  is  the  apogee  of  the 
modified  orbit. 

The  range  angle  traversed  on  the  'nominal  trajectory  is 
from  turnout  to  transfer,  or  On  the  modified  orbit, 

the  traversed  angle  is  from  apogee  to  impact,  or  The 

resulting  range  angle,  is 


“1T 


(33) 


Cost  Function  Parameters 

Tims  of  flight  calculations  are  a"'  in  Chapter  IV.  An 
eccentric  anoaaly  at  transfer  Is  defined  by  using  the  value 

Of 


E2=Cos‘'' 


'gn'CosVp  \ 

\l-ep|Cos  t/g  / 


(34) 
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Since  pre-apogee  transfers  are  considered,  £3  also  lies 
"between  0 and  tt. 


Using  this  result,  the  time  of  flight  on  the  nominal 
trajectory  is  found  from  the  Kepler  time  of  flight  equation 1 

^ ^ 2 ^2 ~ (35) 

where  Tfi2  time  of  flight  from  burnout  to  trajectory 

modification.  Reaction  time,  Tj.,  Is  fotind  from  eq  (25), 

The  transfer  AV  calculations  are  different  from  those  in 
Chapter  IV.  Figure  6 shows  that  the  direction  to  perigee  is 
different  for  the  modified  orbit.  A modified  orbit  perifocal 
reference  frame  is  defined  in  order  to  calculate  velocities 
from  eq  ( 24) , Velocities  on  the  nominal  orbit  are  found  directly 
from  eq  (24)  in  the  nominal  perifocal  reference  frame.  The 
relation  between  the  nominal  and  modified  perifocal  frames  is 


-Cos  1^2 

Sin  1^2 

E 

i 

(35) 

-Sin  1/2 

-Cos  1^2 

_ 

V 

qm 

1-  J 

as  found  from  the  true  anomaly  at  transfer,  i'gt 
geometry  of  the  reference  frames.  From  this  relation,  the 
modified  orbit  velocities,  Vpjg  and  V^j^j  can  be  written  in 
nominal  perifocal  frame  ns  Vpt  and  Vq»  for  direct  comparison 
with  the  velocities  of  the  nominal  trajectory. 
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Froa  eq  { 24)  the  velocity  components  on  the  nominal 
orbit  at  transfer  are 


- Sin 

V^=  - 


1/2 


P [r2(1’enCosP'2?l 


V - Cost-2 
P [^2^^>e„CosU2)f^^ 


(37) 


(38) 


where  Vp  and  are  t .e  velocities  in  the  nominal  perifocal 
fraae. 

Because  transfer  to  the  modified  orbit  occurs  at  apogee, 
velocity  in  the  aodlfied  ? direction  is  zero^  Vpj^  = 0,  But 


V, 


1-e 

m 


qm 


(39) 


as  before.  The  radius  at  transfer,  T2»  apogee  radius 

of  the  new  orbit. 

When  the  modified  orb! t velocity  is  transformed  to  the 
nominal  perifocal  frame  by  eq  (36),  the  transfer  A V is  found 
for  ^ch  nominal  perifocal  direction  i 


^p'  ” ^2  ~ ^p 


(40) 
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By  combining  the  two  differences,  the  total  AV  is  foiindi 


(AVp2.AV^2j1/2 


(42) 


Now,  since  AV  and  Ty  are  defined,  a cost  may  be  found 
for  each  nominal  trajectory  and  its  corresponding  solutions, 

Wethod  of  Solution 

For  the  pre-apogee  to  apogee  transfer,  T2  may  vary 
between  the  burnout  and  apogee  radii  of  each  nominal 
trajectory.  An  infinite  number  of  modified  orbits  which 
hit  a required  final  range  angle  exist  for  a single  nominal 
trajectory.  The  solution  of  this  transfer  problem  must 
find  the  optimal  transfer  for  a given  nominal  trajectory. 

The  method  of  finding  this  optimal  transfer  uses  two 
steps.  First,  the  values  of  rg  and  are  fixed,  and  a 
trajectory  which  connects  the  transfer  point  and  the  target 
is  computed.  For  a constant  a series  of  solutions  is 
found  by  allowing  T2  to  vary  between  burnout  and  apogee. 

A set  of  cost  values  is  then  calculated.  These  cost  values 
are  sorted  to  find  the  lowest  value  of  cost  and  the  value 
of  T2  for  which  It  occurs.  This  solution  then  represents 
the  region  of  a minimum  cost  pre-apogee  to  apogee  transfer 
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from  a given  nominal  trajectory.  Likewise , vvther  transfers 
are  found  by  calculating  solutions  for  different  nominal 
range  angles.  The  res'ilt  of  this  computation  is  a set  of 
near-optimum  pre-apcgA  to  apogee  transfers  which  vary  with 
nominal  trajectory.  Fi' >gram  3 of  Appendix  C is  written  to 
accomplish  this  first  step  of  the  solution. 

The  region  of  a minimum  cost  transfer  is  predicted  from 
the  first  step.  Now,  the  actual  minimum  ccst  trans/'er  for 
a single  nominal  trajectory  is  found  by  solving  tne  first 
order  necessary  conditions  of  a minimum.  The  initial  estimates 
for  the  Iterative  solution  are  obtained  from  the  first  step. 

As  before,  the  augmented  cost  fimctlon,  eq  (23),  is 
■formed.  The  variables  in  this  problem  are  em,  r2»  and  X.  The 
first  order  necessary  conditions  are 


— =C  = AV  .WTr-r^  A — 

‘ 60  ' 60 

m m m 


(4  3) 


^ = 0 = AV  'WT  ‘ 
6 p2  ^ O p2 


O T2 


(44) 


M = 0 A"  - 
dX 


(45) 


The  partial  derivatives  of  J are  found  as  discussed  in 
Chapter  IV,  Equations  (43),  (44),  and  (45)  are  solved  by 
an  itei-ation  on  e^j,  r2,  and  X.  In  this  way,  the  optimum 


pre-apogee  to  apogee  transfer  for  a given  nominal  range 
angle  is  found.  A series  of  solutions,  where  is  allowec 
to  vary,  then  produces  the  cost  function.  Program  4 of 
Appendix  C is  written  to  complete  this  second  step. 

Results  and  Analysis 

For  a target  range  of  100®,  the  minimum  cost  solutions 
were  found  with  nominal  range  angles  varying  between  70° 
and  180°.  A weighting  factor  of  .15  was  employed  in  uhese 
calculations.  Results  are  illustrated  in  figs.  7 through  10. 

As  shown  by  fig,  7,  the  cost  function,  in  this  case, 
decreases  up  to  the  maximum  nominal  range  boundary.  Figure 
8 shows  that  up  to  a certain  nominal,  range  angle  (approximately 
138°),  the  optimum  transfer  point  is  the  nominal  trajectory 
apogee.  Beyond  this  range  angle,  the  optimum  transfer  point 
in  the  nominal  trajectory  moves  toward  burnout. 

The  change  in  the  cost  function,  as  compared  to  the  first 
case,  is  due  mostly  to  the  different  beliavior  in  the  transfer 
AV.  Figures  9 and  10  show  that  as  nominal  range  increases, 
the  reaction  time  remains  a generally  decreasing  function,  but 
the  AV  differs  from  the  first  case.  At  the  point  where  apogee 
to  apogee  transfers  become  non-optimal,  the  AV  begins  to 
decrease.  Up  to  this  point,  AV  is  an  increasing  function. 

This  beha'«'ior  accounts  for  the  decreasing  trend  in  the  cost 
function  at  the  larger  nominal  range  angles. 

The  initial  decrease  of  AV  shown  in  fig.  10  is  due  to 
the  fact  that  the  magnitude  of  AV  is  calculated.  The 
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Vjjj-Vjj  Ts,  range  plot  of  the  apogee  to  apogee  transfer  case, 

5)  shows  that  the  AV  for  the  lower  range  angles  Is 
applied  in  the  direction  opposite  to  vehicle  aotioni 
however,  fig.  10  shows  velocity  sagnitude  only. 

The  decrease  in  AV  at  the  larger  nominal  range  angles 
is  explained  by  intuitive  insight  into  this  transfer  problea. 
For  >iU^  a relatively  high  eccentricity  of  the  aodlflwi  orbit 
is  required  for  apogee  to  apogee  transfers.  If  transfer 
is  done  prior  to  the  nominal  orbit  apogee,  a lower  eccentricity, 
i.e.,  a smaller  orbit  change,  le  adeqmte  enough  to  hit  the 
required  final  range,  These  facts  explain  the  behavior  of  this 
transfer  case. 

The  cost  function  found  for  pre-apogee  transfers  at 
higher  nominal  range  angles  points  out  a significant  advantage 
in  fractional  orbital  ballistic  systems  (?X)BS),  The  vehicle 
is  launched  into  a low,  near  circular  orbit.  Before  com- 
pletion of  the  orbit,  a retro  rocket  slows  down  the  vehicle, 
causing  it  to  drop  on  the  target.  Bamaby  (Ref  1«20)  discusses 
this  type  of  system  and  its  tasic  advatage  is  ttet  the  low 
orbit  remains  imdetected  by  ground  radar  until  the  range  to 
the  target  is  approximately  1400  km. 


VI. 


POST-APOGEE  TO  APOGEE  TRANSFERS 


The  thli^  case  considered  Is  the  post-apogee  transfer 
from  a nominal  orbit  to  apogee  on  the  new  orbit.  Once  again, 
a limiting  case  of  this  transfer  is  the  apogee  to  apogee 
transfer  described  in  Chapter  IV. 

This  transfer  problem  is  essentially  the  same  as  th^; 
problem  of  Chapter  V,  vdth  the  exception  that  the  true  anomaly 
at  transfer  is  after  nominal  apogee.  This  transfer  is 
illustrated  in  fig,  11. 

Modified  Orbit  and  Cost  PHinction  Parameters 

All  parameters  for  thl s case  are  the  same  except  those 
of  true  anomaly  and  eccentric  anomaly  at  transfer.  Equations 
(32)  and  (3^)  are  employed,  but  it  is  noted  that  both  5^2 
ami  E2  lie  between  tt  and  2tt, 

Method  of  Solution 

The  method  of  computing  an  optimum  cost  for  each  nominal 
tra;Jectory  is  as  discussed  in  Chapter  V.  Both  Program  3 
Program  4 of  Appendix  B are  used  to  find  optimum  solutions. 

Results  and  Analysis 

The  cost  function  found  for  this  transfer  case  is 
depicted  In  fig,  12,  Figure  I3  shows  that  for  the  post- 
ap'gee  to  apogee  transfer,  the  optimum  transfer  point  occurs 

.•if ter  nominal  apogee  for  the  lovrer  nominal  range  angles. 

As  range  increases,  ^he  optimum  transfer  point  moves  first 
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towa2*d  the  target,  then  hack  to  apogee,  with  apogee  being  the 
optimum  transfer  point  for  large  nominal  ranges. 

Because  this  is  a post-apogee  to  apogee  transfer,  the 
reaction  time  would  be  expected  to  be  lower  than  the  first 
two  cases.  The  reaction  time  is  shorter  than  for  the  other 
two  cases,  but  fig.  1^  shows  that  the  reaction  time  increases 
up  to  the  point  where  apogee  to  apogee  transfers  become 
optimal.  At  this  nominal  range  angle,  the  reaction  time  starts 
decreasing. 

Figure  15  shows  that  the  transfer  AV  initially  decreases, 
up  to  the  apogee  to  apogee  transfer  point.  Beyond  this 
nominal  range  angle,  AV  Increases  with  nominal  range. 

These  results  can  also  be  intuitively  explained.  As 
would  be  expected,  the  post-apogee  to  apogee  transfers 
require  a higher  eccentricity  in  order  to  bit  the  target 
range.  This  Implif-s  thfit  a higher  AV  than  the  pre-apogee 
case  is  needed^  The  magnitude  of  AV  is  higher  as  expected, 
but  this  higher  impulse  requirement  is  offset  by  a decrease 
In  reaction  time.  The  trade-offs  between  each  cost  component 
as  determined  by  the  weighting  factor„  produce  the  point  at 
>diich  transition  from  post-apogee  transfers  to  apogee 
transfers  becomes  optimal.  The  weighting  factor  also 
produces  tne  behovlcr  exhibited  by  the  true  anomaly  at 
transfer  in  fig.  13. 
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VII.  PROBLEM  DEFINITION  - THREE  DIMENSIONAL  TRANSFERS 


Several  closed  form  mid-course  transfers  have  been 
examined  in  the  preceedlng  sections  to  establish  the  basic 
behavior  of  the  mid-course  modification.  To  see  real  world 
behavior,  a more  complex  model  of  both  the  nominal  and 
modified  trajectories  is  needed. 

Nominal  Trajectory  Definition 

As  previously  mentioned,  the  typical  ballistic  missile 
trajectory  is  lofted.  A real  world  simulation  should  use 
this  type  of  nominal  trajectory.  In  addition,  earth  rotation 
must  be  taken  into  consideration.  For  a target  at  45®  latitude, 
the  impact  error  due  to  earth  rotation  is  on  the  order  of 
400  miles  for  a normal  ballistic  missile  trajectory  (Ref  8*2?). 

Earth  oblateness  affects  the  trajectory  by  changing  the 
gravitatior.al  accelerations  from  motion  in  an  ideal  inverse 
square  gravitational  field  to  accelerations  which  are 
dependent  on  latitude.  For  ICBM  ranges,  errors  caused  by 
ignoring  oblateness  effects  are  on  the  order  of  10  miles 
(Ref  8*36).  Earth  oblateness  not  only  affects  vehicle 
accelerations,  but  initial  and  final  position  errors  result 
when  locating  a point  on  a spherical  earth  instead  of  using 
an  oblate  spheroid  model.  Position  errors  are  latitude 
dependent  and  may  vary  to  as  much  as  21.4  km,  as  this  is 
the  measured  equatorial  bulge  of  the  earth  (Ref  2*95). 

Therefore,  it  is  imperative  that  earth  oblateness  be 
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considered  in  a realistic  model 


Finally,  atmospheric  entry  of  the  reentry  vehicle  must 
be  simulated.  Aerodynamic  deceleration  of  the  reentry  vehicle 
directly  affects  the  trajectory  and  must  be  included.  A 
numerical  simulation  which  accounts  for  all  these  real  world 
effects  is  used  in  this  second  pert  of  the  investi^tion. 


Assumptions 

To  generate  a trajectory,  boundary  conditions  at  bum- 
out  and  the  target  location  must  be  established.  The  position 
of  the  missile  with  respect  to  the  launch  site  at  burnout 
locates  the  burnout  point  in  an  earth  centered  inertial  frame. 
Velocity  imparted  to  the  missile  by  the  rocket  motor  is  added 
to  the  launch  site  inertial  velocity  to  determine  missile 
inertial  velocity  at  burnout.  Several  assumptions  are  made 
for  these  initial  boundary  conditions.  The  magnitude  of  a 
radius  vector  from  the  launch  site  to  the  missile  at  burnout 
is  set  at  200  n.  ml.  The  elevation,  and  azimuth  from  north 
of  the  radius  vector  from  the  launch  site  to  the  burnout 
point  are  assumed  to  be  45°  and  20°  respectively.  The 
magnitude  of  the  velocity  imj^rted  to  the  missile  is  assumed 
to  be  25.000  ft/sec,  a typical  ballistic  missile  burnout 
velocity. 

The  location  of  the  launch  site  and  target  on  an  oblate 
earth  are  necessary  to  establish  boundary  point  positions. 
These  locations  are  defined  in  terms  of  lat'tude  and  longitude 
angles.  The  launch  site  and  target  positions  are  needed  for 
an  earth  rotation  model i Inertial  velocity  imparted  to  the 
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missile  by  the  launch  site  varies  with  launch  site  latitude. 
The  movement  of  the  target  in  inertial  space  is  also  latitude 
dependent. 


Oblate  earth  effects  upon  vehicle  motion  are  considered 
by  including  four  gravity  harmonics  (Ref  2i4l9),  These 
harmonics  account  only  for  gravity  anomalies  symmetrically 
distributed  about  the  earth's  spin  axis.  Sectoral  and  tessera  1 
harmonics,  dependent  on  specific  longitudes,  are  not  included 
because  their  effects  vary  for  different  trajectories. 
Trajectory  perturbations  caused  by  the  sun  and  the  moon  are 
small  and  so  are  not  included  (Ref  6:42,  8!45). 

Adding  reentry  to  the  model  requires  several  assumptions. 
An  exponential  atmospheric  model  is  used  for  density  calcu-' 
latlons.  Atmospheric  effects  are  used  only  when  the  altitude 
is  below  approxl^tely  110  km  (Ref  6:44),  Altitude  compu- 
tations assume  a flat  earth  in  the  vicinity  of  the  target 
because  the  range  covered  by  the  reentry  vehicle  after  reentry 
is  sJiBll  (Ref  6:A4),  The  altitude  is  then  found  by  finding 
the  difference  between  magnitudes  of  the  inertial  radius 
vector  to  the  vehicle  and  the  inertial  radius  vector  to  the 
target.  This,  in  effect,  produces  an  oblatt:  atmospheric  model. 

For  drag  calculations,  some  measure  of  vehicle  stream- 
lining is  needed.  A ballistic  coefficient,  relating  vehicle 
imsSf  frontal  area,  and  drag  coefficient  to  the  total  drag 
is  employed.  The  use  of  this  quantity  in  drag  equations  is 
defined  in  Appendix  A.  The  value  of  the  ballistic  coeniclent 
used  is  4500  kg/i3“  corresponding  to  a streamlined  or  heavy 
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reentry  vehicle 


The  fact  that  the  ataosphere  rotates  in  inertial  space 
is  also  considered  in  the  drag  calculations.  The  rotation 
of  the  atmosphere  has  a net  effect  of  changing  the  velocity 
of  the  vehicle  with  respect  to  the  air.  Atmospheric  drag 
acts  as  a force  opposite  to  the  direction  of  velocity. 

Direct  entry  into  the  atmosphere  is  assumed,  with  no  skipping 
reentry  considered.  Additionally,  no  lift  due  to  vehicle 
attitude  is  Included, 

Finally,  a time  of  powered  flight,  between  missile 
launch  and  burnout,  is  needed  to  determine  how  far  the 
launch  site  and  ^■arget  have  moved  during  the  powered  phase. 

Time  of  powered  flight  is  assumed  to  be  5 minutes. 

The  assumptions  made  for  the  nominal  trajectory  derivation 
'vfe  summarized  in  Table  I. 


T.iiBLE  I 


LOFTED  NOMINAL  TRAJECTOKV  ASSUMPTIO.NS 


Parameter 


Burnout  Position  Vector* 
Distance 
Elevation 
Azinnath 


Burnout  Velocity  Magnitude* 


Reentry  Altitude 


Vehicle  Ballistic  Coefficient  i 


Time  of  Powered  Flight 


* With  respect  to  the  latinch  site. 


Assumed  Value 


200  n mi 
20° 


25,000  ft/sec 


110  km 


4500  kg/ 


5 min 
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The  use  of  the  proceeding  assumptions  in  defining 
boundary  points  and  the  equations  of  motion  for  the 
numerical  simulation  of  the  trajectory  is  discussed  further 
in  Appendix  A, 

Method  of  Solution 

Because  the  target  is  moving,  three  tine  varying  Inertial 
coordinates  determine  final  conditions.  To  satisfy  the  three 
terminal  conditions  of  the  target,  three  Initial  conditions 
at  burnout  are  adjusted.  The  velocity  magnitude  is  fixed, 
so  only  the  direction  may  be  varied  at  burnout.  The  burnout 
point  is  fixed,  and  only  one  other  parameter,  time  of  flight 
to  impact,  remains.  The  three  variables,  azimuth  of  the 
velocity  vector.  Its  elevation,  and  time  of  flight,  can  be 
varied  at  burnout  so  that  the  vehicle  Impacts  upon  the  target 
at  the  final  time. 

The  three  variables  are  used  in  Powell’s  algoritrim  and 
the  three  functions  are  formed  by  a three-dimensional  "hit" 
condition  at  final  time: 


=0  = X,-X^ 

(46) 

F2=0=Yf  -Y^ 

(47) 

F3  = 0 = Zf-Z. 

(48) 

where  X,  Y,  and  Z are  the  inertial  coordinates 
f refers  to  the  vehicle  at  fiml  time 
t refers  to  the  target  at  final  time. 
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To  Insure  convergence  to  the  lofted  trajectory,  the 
Initial  estiaate  of  the  velocity  vector  elevation  is  set 
at  60°,  Equations  (^),  (4?),  and  (48)  are  then  used  in 
Powell's  routine  where  the  boundary  conditjcns  derived  in 
Appendix  A determine  trajectory  end  points.  The  vehicle 
behavior  in  inertial  space  is  modeled  by  the  nonlinear 
equations  of  motion,  also  found  in  Appendix  A.  The  equations 
of  motion  are  integrated  numerically,  using  a Runge-Kutta 
and  four  point  predictor-corrector  library  routine.  The 
trajectory  Initial  conditions  are  varied  by  Powell's  routine 
so  that  eqs  (46),  (47),  and  (48)  are  satisfied.  The 
solution  of  a nominal  trajectory  is  incorporated  in  Program 
5 of  Appendix  C. 
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(53) 


— — -0  = Av  ♦X,- — i— *Xp  — ^ — • X-5— J— 

dAV^  ^ ’dAV^  ^ 6AV..  -^6AV, 


61^; 

_5Y^ 

dTf 

,s  1 \ 

(54) 


=0  = X,  - X. 
6Xi  f ^ 


(55) 


6J 

dX 


= 0 = Y^-Y^ 


(56) 


^ =0=  Z,-2, 

6a-^  • ^ 


(57) 


These  equations  may  now  be  solved. 

Method  of  Solution 

The  solution  to  this  problem  provides  conditions  at 
a point  of  modification  in  the  nominal  trajectory.  The 
Initial  conditions  of  the  modified  trajectory  define  r ' 
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trajectory  which  hits  the  target  and  also  represents  a 
minlamm  cost  transfer  from  the  point  of  trajectory  change. 

The  partial  derivatives  of  J with  respect  to 
AVg  and  T^,  are  needed.  These  are  found  numerically  by 
calculating  first  differences. 

Since  there  are  seven  equations  end  seven  unknowns, 
Powell’s  alg'-rithm  can  be  employed  to  satisfy  the  first  order 
necessary  conditions  of  a minimum.  The  Initial  values  of 
AV^,  AVyf  AV^  and  T^,  can  be  roughly  estimated  by  generating 
two  nominal  trajectories.  The  first  is  between  the  launch, 
site  and  pseudo  target.  The  second  is  from  the  launch  si*;e 
to  the  real  target.  A time  of  flight  along  the  nominal 
trajectory  defines  the  point  of  transfer.  The  initial 
estimates  for  AV^,  AVy,  and  AV2  are  found  by  calculating  the 
difference  of  the  velocity  components  of  the  two  trajectories 
at  the  time  of  transfer.  The  estimate  of  T^.  is  the  difference 
between  the  total  time  on  the  second  trajectory  and  the 
time  at  modification.  The  Lagrange  multipliers  are  initially 
set  to  zero.  The  method  of  finding  initial  estimates  is 
presented  in  Table  II,  These  rough  initial  estimates  are 
within  the  area  of  convergence  for  the  algorithm  used. 

The  same  equations  of  motion  used  for  the  nominal 
trajectory  generation  are  employed  for  the  modified  orbit 
calculation.  The  equations  of  motion  are  integrated 
numerically  to  solve  eqs  (51)  through  (57).  This  method 
t^as  programmed  and  .is  incorporated  in  Pi-ogram  5 of 
Appendix  C. 
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TABLE  II 


INITIAL  ESTIMATES  - NONPLANAR  TRAJECTORY  MODIFICATION 


Variable 

j Initial  Estimate  at  Time  of  Modification 

AVi 

Vxm  “ ^xn 

modified  orbit  Vx  - nominal  orbit  Vx 

AVy 

Vym  - Vyn 

modified  orbit  Vy  - nominal  orbit  Vy 

AVz 

^zm  “ 

modified  orbit  Vg  - nominal  orbit 

Tr 

^ffm“Tmod 

modified  orbit  free  flight  time 
- time  at  transfer 

0 

>■ 

0 

n 

^3  1 

0 

Results  i 

To  test  this  method,  a typical  ICBM  range  to  real  and  I 

pseudo  targets  Is  selected.  The  launch  site  is  located  at  j 

3?. 5°  N»  latitude  and  125^  W.  longitude.  The  coordinates  of  | 

the  pseudo  and  real  targets  are  N,  latitude,  3.3°  W,  | 

I 

longitude  and  55.5°  N.  latitude,  5.^°  E.  longitude,  respectively,  | 

I 

The  distance  between  pseudo  and  real  targets  is  approxiriiately  i 

350  n mi.  Modification  takes  place  approximately  32  minutes  I 

into  the  nominal  trajectory  62  minutes  of  flight,  f 


The  initial  state  estimates  are  scaled  to  lie  between 
-1  and  1 for  better  performance  of  the  algorithm  (Ref  7il20), 
In  addition,  the  matrix  of  dP/dX,  the  Jacobian,  is  examined 


to  check  that  the  scaling  of  both  F ard  X do  not  produce  an 
ill-conditioned  matrix.  For  this  reason,  the  hit  conditions, 
eqs  (55).  (56),  and  (57)  are  also  scaled. 
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After  scaling  both  F and  X,  the  algorithm  readily 
converged  x:o  a point  where  the  sum  of  the  squares  of  F 
Is  about  three.  At  this  point,  the  algorithm  continues  to 
converge,  but  very  slowly.  In  examining  the  initial  conditions 
being  adjusted,  the  velocity  at  modification  is  changing  by 
about  .0001  ft/sec.  Time  of  flight  is  changing  by  .001 
seconds.  Further  reduction  in  the  sum  of  the  squares  error 
produces  no  practical  results  when  noting  the  initial  condition 
magnitude  changes.  The  fact  that  the  three  dimensional  hit 
eqmtion  is  scaled  not  only  produces  better  convergence,  but 
also  insures  that  the  real  target  is  reached  with  reasonable 
accuracy  when  the  algorithm  is  stopped.  For  the  trajectory 
under  consideration,  the  miss  distances  as  determined  by 
eqs  (55)»  (56) » and  (57)  in  each  inertial  direction  when  the 
algorithm  stopped  are  -497,3  ft,  -510  ft,  and  415,4  ft  in 
the  X,  Y,  and  2 inertial  directions  respectively.  This  is 
translated  into  a spherical  miss  at  the  target  of  ,136  n mi. 

The  magnitude  of  the  transfer  AV  is  I255  ft/sec  and  the 
reaction  time  is  2055.8  seconds.  The  corresponding  cost  for 
this  mid-course  modification  is  ,488  . These  results  are 
summarized  in  Table  III. 

An  Important  consideration  in  selecting  a trajectory  is 
its  sensitivity  to  initial  errors.  Velocity  errors  in  each 
inertial  direction  were  added,  one  at  a time,  and  the 
equations  of  motion  were  integrated  foreward  to  final  time 
in  order  to  determine  : inal  miss  distance  due  to  initial 
velocity  errors.  Table  IV  presents  the  results  of  these 
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calculations.  Table  IV  shows  that  for  this  test  case, 
realistic  velocity  errors  at  the  point  of  orbit  modification 
produce  only  small  final  miss  distances. 

TABLE  III 

NONPLANAR  TRAJECTORY  MODIFICATION  RESULTS 


1 1 

Parameter 

Value 

Launch  Site  Location 

37.5°  N.lat.  125°  W.long. 

Pseudo  Target  Location 

54°  N.lat.  3.3°  W.long. 

Real  Target  Location 

55.5°  N.lat.  5.4°  E.long. 

Velocity  Vector  at  Burnout* 

Azimuth 

19° 

Elevation 

49° 

Time  of  Flight  at  Transfer 

32  min 

{AVI 

1255  ft/sec 

446.3  ft/sec 

1121.8  ft/sec 

AVz 

342.8  ft/sec 

Tr 

2055.8  sec 

J 

.488 

* Burnout  condition  assumptions  are  found  in  Table  I.  These 
values  are  the  result  of  the  iteration. 
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TABLE  IV 

VELOCITY  SENSITIVITY  - NONPLANAR  TRAJECTORY  MODIFICATION 
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Conclusions 

This  study  has  analyzed  the  single  Impulse  mid-course 
ballistic  trajectory  modification.  Methods  for  examining 
both  coplanar  and  three-dimensional  transfers  were  presented. 

Application  of  the  methods  shown  ro  real  systems  is  dependent 
upon  knowing  trajectory  parameters  which  have  been  assumed 
for  illustration  of  the  methods. 

Coplanar  transfers  for  elliptic  orbits  vrlth  no  atmospheric 
reentry  or  earth  rotation  were  first  examined  to  define  the 
general  behavior  of  the  transfer  problem.  The  nominal  orbit 
parameter  of  interest  is  free  flight  range  angle.  Because  a 
maximum  range  nominal  trajectory  is  used,  this  range  angle, 

i 

along  with  the  altitude  of  the  burnout  point,  defines  all  " 

nominal  orbit  parameters.  The  optimum  transfers  for  the 
three  cases  investigated  are  plotted  against  nominal  range 
angle  in  fig.  16. 

As  fig,  16  indicates,  the  pre-apogee  transfer  to  a near 
circular  orbit  produces  the  lowest  cost  for  this  formulation. 

It  was  assumed  that  a booster  is  capable  of  reaching  any 
range  angle  up  to  180°.  If  rhls  is  not  the  case,  the  best 
type  of  transfer  to  use  can  be  found  from  fig,  l6  by  noting 
the  maximum  attainable  nominal  range  angle. 

Optimum  solutions  for  the  coplanar  orbit  transfers  were 
found  by  solving  first  order  necessary  conditions  of  a minimum 
in  the  cost  function.  This  procedure  dees  not  assure  that  a 


61 


111 


Conparioon  - Coplanor  Transfers 


minimum  is  actually  calculated.  This  fact  was  considered 
and  for  this  reason,  the  minimum  cost  solutions  of  Program 
3 were  used  as  a starting  point  to  find  the  minimum.  By 
starting  the  algorithm  in  the  area  where  a series  of  solutions 
show  the  minimum  to  exist,  convergence  to  the  minimum  point 
is  insured. 

The  mid-course  transfer  was  then  extended  from  the 
coplanar  transfer  to  the  case  where  the  nominal  and  modified 
orbits  are  nonplanar.  A model  which  uses  a rotating,  oblate 
earth,  atmospheric  reentry,  and  a rotating,  oblate,  exponential 
atmosphere  is  employed, 

A method  of  finding  the  conditions  necessary  for  an 
optimum  trajectory  modification  was  then  developed.  The 
region  of  a minimum  cost  solution  was  found  by  terminating 
the  algorithm  before  the  first  order  necessary  conditions 
were  satisfied  exactly.  For  the  example  considered,  the 
routine  was  teraJinated  when  velocity  change  at  the  point  of 
trajectory  modification  changed  only  slightly,  beyond 
reasonable  accuracy  of  control  in  a realistic  system.  Final 
miss  distance  at  the  target  still  remained  within  an 
acceptable  distance  when  the  algorithm  stopped.  Practical 
considerations  of  computer  calculating  tine  also  entered  into 
the  decision  of  ■*  minating  the  routine.  Generation  of  the 
one  trajectory  modification  found  for  this  case  took 
approximately  ISOO  seconds  of  computer  time  from  the  rough 
initial  estimates.  Convergence  at  the  termination  point 
was  slow. 
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The  sensitivity  of  the  modified  trajectory  to  initial 
velocity  errors  was  examined  by  perturbing  the  transfer 
velocity  impulse.  It  was  found  that  for  the  modified 
trajectory  calculated,  realistic  velocity  errors  produce 
no  significant  final  miss  at  the  target. 

It  should  be  noted  that  several  ass*amptions  of  launch 
site  and  target  locations,  burnout  states,  and  vehicle 
parameters  have  been  made.  These  are  for  illustration 
purposes  and  application  of  the  method  to  a specific  problem 
will  change  their  values. 

In  addition,  the  cost  functjon  weighting  factor  was 
chosen  to  produce  required  AV  in  a realistic  range  of  values 
The  formulation  of  the  problem  with  a factor  which  can  be 
changed  to  accommodate  trade-offs  between  AV  and  reaction 
time  allows  a trajectory  designer  to  apply  this  algorithm 
to  different  AV  ranges. 

Recommendations 

For  the  elliptical  orb^t  coplanar  transfer,  only 
transfer  to  a modified  orbit  apogee  was  considered.  This 
could  be  extended  to  analyze  transfer  from  a nominal  orbit 
to  a coaxial  elliptic  orbit,  and  then  transfer  to  a 
generalized  orientation  elliptical  crbit.  Because  direction 
of  orbit  change  from  nominal  to  modified  trajectory  is  quite 
different  when  the  target  and  nominal  ranges  vary  greatly, 
it  seems  that  these  other  two  types  of  transfers  could  do 
better  in  terms  of  cost  for  some  nominal  trajectories. 


The  effects  of  changing  the  cost  fxmctlon  weighting 
factor  should  be  considered.  For  values  near  the  example 
weighting  factor  used,  the  cost  function  curves  appear  to 
shift  as  new  optimums  are  defined.  However,  for  different 
values  of  W,  the  cost  function  appears  to  change  appreciably. 
This  behavior  was  not  examined  in  detail,  and  further  study 
is  required. 

Finally,  the  performance  of  an  actual  system  could  be 
estimated  using  the  second  part  of  the  investigation. 

Actual  system  parameters,  which  have  been  estimated  for 
this  study  would  be  needed.  A set  of  optimal  nonplanar 
modified  trajectories  can  be  generated  and  the  one  that 
performs  best,  i.e.,  in  error  sensitivity  and  cost,  can  be 
chosen. 
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frame  through  a geodetic  latitude  angle.  L.  and  an  equatorial 
rotation  angle.  B . similar  to  a local  sideral  time.  The 
roierence  frames  and  their  relations  are  shown  in  fig,  17, 


launch  Site  and  Target  Locations 

The  launch  site  and  target  locations  in  the  inertial 
reference  frame  are  found  through  use  of  a geodetic  latitude, 
an  equatorial  rotation  angle,  and  an  oblate  earth  model, 

where  both  target  and  launch  site  are  assumed  to  be  at  sea 
level. 


La^?lch  Site  Inertial  Position 

The  oblate  earth  model  of  Bate  (Ref  2xiQ)  provides  1 
and  k locations  of  the  launch  site  at  time  of  launch i 


V 1 o. 

(1-e|sin2g1/2 

v|  c ■ aeSinL|(1-ep2) 

(l-e|siri2L|_)'i/2 


(58) 


(59) 


where  XLS  is  the  initial  1 coordinate 
is  the  initial  k coordinate 
®e  is  the  eccentricity  of  the  oblate  earth  (,08l8l) 
Ll  is  the  launch  site  geodetic  latitude 
ag  is  the  equatorial  radius  of  the  earth. 
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If  the  r and  3 coordinates  at  any  time  after  launch  are 
desired,  the  rotation  angle,  is  employed.  For  the  launch 
point 

0|_=WeT  (60) 


where  Wg  is  the  earth’s  rotation  rate 
T is  the  time  of  missile  flight. 

The  k component  of  the  lavinch  site  does  not  change. 
Therefore  the  Inertial  location  of  the  launch  point  at 
any  time  is 


”XLS  CosO^ 

= 

XLS  Sin  (9, 

ZLS 

LS^l 

Rr 


(61) 


LS_I 

where  Rj  Is  the  la\mch  site  location  with  respect  to 
inertial  frame,  written  in  the  inertial  frame. 


Target  Inertial  Position 

Similarly,  the  target  location  in  the  inertial  reference 
frame  can  be  found.  Define  XTS  and  ZTS  such  that 


XTS 


a^Cos  Lj 
( 1 - e^Sin"^ Ly 


(62) 
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(63) 


ZTS  - n Lj(  1 ~ ) 

(1-e|sin2LT.)'’/2 

where  Lp  is  the  target  geodetic  latitude. 

The  equatorial  rotation  angle,  is  found  once  knowing 
the  time  and  the  Initial  longitude  difference  between  target 
aivi  launch  slte^  N. 

0-p=cOgT*N  (64) 


The  angle  N is  measured  eastward  from  the  launch  point. 
The  inertial  position  of  the  target  at  an^’’  time  is 
therefore 


XTS 

XTS 


Cos  Oj 
Sin 


ZTS 


(65) 


T_I 

where  Hj  is  the  target  location  with  respect  to  the  inertial 
frame,  written  in  the  inertial  frame. 


Burnout  Point  Inertial  Location 

To  find  the  bvimout  point  in  inertial  s|ace,  it  Is  first 

necessary  to  define  & vector  from  the  launch  y>oint  to  the 
burnout  point  at  time  of  burnout.  This  vector,  written  in  the 
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launch  site  topooentric  frame  can  he  defined  If  a distance, 
azlcmth  and  elevation  (corresponding  to  a radar  measurement) 
ere  Known.  The  vector  is  then  written  In  the  launch  site 
topooentric  frame  at  burnout i 


- Cos  iSpj  Cos^^p  R[jq 
Sin^Sp  Cos4>p  R(30 

Sin<?>R 


(66) 


where  Is  the  azimuth  angle  from  true  north  of  the  burnout 
vector 

Is  the  elevation  angle  of  the  vector 
B^o  Is  the  magnitude  of  the  vector. 

The  topooentric  fiame  at  the  Isxinch  site  is  related  to 
the  defined  Inertial  reference  frame  by  two  Euler  angle 
rotations  through  a latitude  angle,  Lj^ , and  a rotation  angle, 

. The  resulting  transformation  )natrix  from  l)he  topooentric 
to  the  Inertial  frame,  [Lj^] , Is 


I 


Sin  L^Cos^^ 
Sin  L^Sin 


-Cos  L, 


-S  i n 9^  Cos  L ^Co  s (9^ 
Cos  9^  CosL^Sin^L 
0 SinL^ 


(67) 
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i^ere  is  v..efined  by  eq  (60),  It  Is  noted  that  a time 
from  latmoh  to  burnout  is  needed  to  determine 

Using  this  transformation,  along  with  the  vectors  which 
describe  the  burnout  point  with  respect  to  inertia,  eqs  (66) 
and  (6l),  the  burnout  point  is  located  in  the  inertial 
reference  frame i 


bo-I 


= l4t] 


bO-L  L5-I 


bo_I 

i^ere  Rj  is  the  vector  from  the  earth's  center  to  the 
burnout  point,  written  in  the  inertial  frame.  This  vector 
gives  the  initial  position  of  the  trajectory. 


Initial  Inertial 


THe  initial  inertial  velocity  is  also  required  to 
complete  the  set  of  initial  conditions  of  the  trajectory. 
This  velocity  is  found  in  a method  similar  to  t»»t  for  the 
position  at  burnout. 

If  a velocity  azimuth  angle,  elevation  angle, 
and  a magnitude  are  known,  the  velocity  vector  in  the 
topocentric  frame  can  be  written  as 


bO-LS 

V 

T 


'Cos  '^bo 

Sin  /9yCoscJ>v 

S in  VjjQ 


It  Is  noted  that  the  angles  and  are  not  related 
to  the  angles  and  which  describe  the  burnout  point  with 
respect  to  the  launch  site.  However,  these  angles  should  be 
sosei^t  similar  as  seen  by  a physical  examples  A missile 
would  not  reach  a bumout  point  west  of  a launch  site  with 
an  easterly  velocity. 

At  this  point,  it  is  necessary  to  note  that  an  initial 
inertial  velocity  is  Imparted  to  the  missile  by  the  launch 
site.  This  inertial  velocity  is  due  to  the  earth’s  rotation 
and  varies  with  the  latitude  of  the  launch  site.  The  velocity 
Is  directed  along  the  easterly  axis  of  the  local  topocentrlc 
frames 


LS-I 

V = VgCos  e 


(70) 


where  V_  is  the  eouatorlal  tangential  velocity  (1524  ft/sec) 
eg  is  a unit  vector  in  the  easterly  direction. 

By  combining  eqs  {69)  and  (70)  and  employing  eq  (67), 
the  initial  velocity  can  be  written  in  the  Inertial  frames 


bo -I 

V, 


/bO-LS  LS-I 


T 


V 


(71) 


The  Equations  of  Motion 


The  differential  equations  of  motion  in  the  inertial 
reference  frame  are  required  to  complete  the  general  trajectory 


definition.  Two  contributions  to  the  acceleration  experienced 
by  the  vehicle  are  examined.  They  are  the  gravitational 
accelerations  and  drag. 


Gravitational  Accelerat 1 ons 

A gravity  model  which  employs  four  gravitational 
harmonics  is  used,  (Raf  2i421)j 


f 
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2\R 
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R / 


Z 
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231 


(72) 


Y •• 


(73) 


75 


« -1.6  X 10“^ 

« ,15  X 10"^ 

This  gravitational  oodel  acoounts  only  for  gravity 
anoaalies  symaetrically  distributed  about  the  inertial  1c 
axis.  Tesseral  and  sectorial  harmonics,  along  with  the 
pertitrbative  effects  of  other  h^venly  bodies  such  as  the 
sun  and  moon  are  Ignored  since  only  small  errors  result 
{Ref  6»i^2,  8i^5). 

The  acceleration  caused  by  drag  is  also  considered  in 
the  equations  of  motion.  It  is  assiuaed  that  burnout  occurs 
above  the  earth's  atmosphere,  so  only  during  reentry  Is  drag 
considered.  The  atmospheric  entry  angle  ias  not  calculated, 
and  It  is  assumed  that  the  vehicle  enters  the  atmosphfere 
directly,  %rlthout  skipping. 

The  acceleration  caused  by  drag  (Ref  5*39) » ®ay  then 
be  written  as 


V_I 

d 


PCqA 


I 


2m 


(75) 


V I 

where  Is  the  drag  acceleration  in  the  inertial  fwme 
P is  the  local  atmospheric  density 
bp  i 1 thu  v^jhiole  irag  coefficient 
A is  the  projected  frontal  a.r^ 

B is  the  ^enicle  mass 
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V is  the  asgaitude  of  the  velocity  vector 

Vj  is  the  inertial  velocity  vector. 

This  acceleration  acts  directly  along  the  flight  jsath 
and  opposite  to  the  velocity  vector. 

The  coefficient  CjvA/a  is  dependent  upon  the  vehicle, 
and  its  Inverse  is  call«i  the  ballistic  coefficient,  BCi 

BC=7t^  (76) 

CqA 

The  ballistic  coefficient  Is  a seasure  of  vehicle  streaa> 
lining.  High  values  of  the  ballistic  coefficient  indicate 
a streaallned  vehicle. 

An  atnospheric  aodel  is  needed  to  find  the  local 
atmospheric  density  required  in  eq  (75)»  An  exponential 
density  model  {Ref  5*38)  is  used,  and  atmospheric  effects 
are  oensidered  negligible  above  an  altitude  of  110  xa.  for 
c<»patatlonal  efficiency^ 

The  fact  that  the  mrth's  atmosphere  rotates  in  inertial 
space  can  alao  be  considered  in  the  di«g  model  (Ref  ZikZk), 
This  phenomena  affecta  the  drag  by  changing  the  velocity 
relative  to  the  atmosphere.  If  the  rotating  ataosphara  la 
include,  the  inertial  directicm  drag  expressions  are 

•*  Pg 
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'^d  = 


exp(-h/23999.3)(V -U)„X)V 


(78) 


Zh=  exp( -h/ 23999.3 )V,V 

d 2BC  2 


(79) 


wh«r«  X^,  and  are  the  inertial  drag  accelerations 
( f t/sec^ ) 

pQ  ia  aea.  level  density  (slug/ft^) 

h is  the  local  altitude  (feet) 

'7  , Y , and  V are  the  inertial  velocity  coaponents 
X • z 

(ft/sec) 

V is  the  velocity  vector  aagnitude. 

The  local  altitiaie,  h,  is  found  by  assttain.i  a flat 
earth  in  the  violnlty  of  the  target.  Because  the  range 
covered  during  reentry  is  relatively  short,  this  is  a 
valid  assuaption  (Hef  6t44}.  Altitude  is  then 


h = 


V ^ 5 1 -r  T ! 

▼ * I If  — — A I 


R. 


I R,i 


(80) 


where  Hj  is  the  i»gnltude  of  the  radius  vector  to  the 
vehicle  with  respect  to  the  inertial  fraae 

T_I 

is  the  aagnitude  of  the  target  zadius  vector  in 
the  inertial  frase.  This  i^gnitude  is  found  froa 
the  oblate  earth  aodel  described  by  eq  (65). 
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By  coablning  eqe  (?2),  (73).  (74),  (77),  (78),  and 
(79).  the  differential  equations  of  notion  result i 


• • 

>^g-Xd 
• • • » 

(81 ) 

Y = 

Y * V 

‘d 

• • • • 

(82) 

Z = 

(83) 

These  differential  equatlcma  of  notion  can  then  be 
Integrated  nunerically  to  find  Inertial  positions  and 
velocities. 


Appendix  B 


THE  NONLINEAfi  BWTIC^  SOLVES 

The  results  of  this  study  are  dependent  upon  solving 
a set  of  highly  nonlinear  botindary  value  prohless.  The 
problems  say  be  stated  as  n functions  of  n unknowns.  In 
order  to  solve  this  set  of  defined  problems,  an  algorithm 
for  solving  systems  of  nonlinear  algebraic  equations  devel- 
oped by  M,  J.  D.  Powell  (Hef  7j115)  Is  employed.  This 
algorithm,  named  subroutine  HSOlA  In  the  computer  programs, 
solves  the  system  of  equations  F(X),  where  ? and  X are  of  the 
same  dimension.  The  Fowcll  algorithm  uses  a modified  Newton 
iteration  where  the  steepest  descent  of  ?(X)  is  also  considered. 
Excellent  results  were  obtained  in  using  this  routine. 

Care  must  be  taken,  however,  to  Insure  that  correct  scaling 
of  both  f and  X is  accomplished.  With  correct  scaling,  the 
algortha  usually  converges  for  even  poor  initial  estimates. 

Both  bo-ondary  value  and  minimization  problems  were  solved 
through  use  of  the  routine. 
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Appendix  C 


COMPUTER  PROGRAM  LISTINGS 

This  section  presents  the  computer  programs  used  in 
this  study.  The  programs  are  written  in  Fortran  Extended 
for  use  on  the  CDC  6600  series  computer. 
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CALCULATE  THE  PARTIALS  NRT  ECCENTRICITY  OF  MODIFIED  ORBIT 
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HAS  THE  VEHICLE  ENTERED  THE  ATMOSPHERE? 
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